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ABSTRACT 



We place firm upper limits on the global accretion history of massive black holes at z > 5 from the recently measured unresolved 
fraction of the cosmic X-ray background. The maximum allowed unresolved intensity observed at 1.5 keV implies a maximum 
accreted-mass density onto massive black holes of p acc < 1.4 x 10 4 M Q Mpc~ 3 for z > 5. Considering the contribution of lower-z 
AGNs, the value reduces to p acc < 0.66 x 10 4 M Q Mpc~ 3 . The tension between the need for the efficient and rapid accretion required 
by the observation of massive black holes already in place at z > 7 and the strict upper limit on the accreted mass derived from the 
X-ray background may indicate that black holes are rare in high redshift galaxies, or that accretion is efficient only for black holes 
hosted by rare galaxies. 
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1. Introduction 

While there is ample evidence that some supermassive black 
holes with masses exceeding 10 9 M formed as early as z > 6 
(with a redshift record of z = 7.1 reported by Mortlock et al. 
2011; see also, e.g., Fan et al. 2001, Barfh et al. 2003, Willott 
et al. 2009, Jiang et al. 2009), there is currently little or no con- 
straint on the evolution of the supermassive black-hole popula- 
tion, as a whole, at the same redshifts. We have only been able 
to probe the most exceptional quasars, which are powered by the 
most-massive black holes. 

A powerful tool capable of globally constraining the nature 
of the high redshift massive black-hole (MBH) population, at 
least of its active fraction, that is manifested as active galactic 
nuclei (AGNs), is the measure of the unresolved cosmic X-ray 
background (CXRB) (Dijkstra et al. 2004, Salvaterra et al. 2005, 
Salvaterra et al. 2007, McQuinn 2012). 

Chandra deep observations have succeeded in resolving al- 
most the entire (80-90%) CXRB over its whole X-ray bandwidth 
(0.5-8 keV). The resolved fraction is almost 100% at low ener- 
gies, but decreases slightly, down to ~ 85%, at higher energies 
(see fig. 8 in Moretti et al. 2012, see also Moretti et al. 2003, 
Worsley et al. 2005). Cosmic X-ray background sources have 
been found to be mostly AGNs with some contribution at soft 
energies (< 2 keV) from galaxy clusters and starburst galaxies 
(Xue et al. 2011, Lehmer et al. 2012). Most of the CXRB sig- 
nal comes from sources located at z < 2, with only ~ 1% being 
produced at z > 4 (Xue et al. 201 1). While in the hard band (2- 
10 keV) the residual unresolved fraction is commonly believed 
to be entirely due to the integrated emission of undetected point 
sources, in the softer band (0.5-2 keV) most of the diffuse emis- 
sion is due to thermal radiation from the Galaxy and the local 
hot bubble (Kuntz & Snowden 2000). A direct assessment of the 
unresolved fraction of the CXRB was performed by Hickox & 
Marckevitch (2007) using Chandra deep field data. They found a 



small but statistically significant diffuse emission in the 1-2 keV 
band, but an emission consistent with zero at higher energies. 
However, the high Chandra instrument background, -25 times 
higher than the unresolved CXRB, makes this measure highly 
uncertain. 

Moretti et al. (2012) exploited the very low (compared to 
Chandra) instrument background of the Swift XRT to measure 
the unresolved spectrum with the highest accuracy available to- 
day. This spectroscopic measure allowed the unresolved CXRB 
to be accurately probed with a much higher energy resolution. 
In particular, the constraint on the 1.5-2 keV band is a factor of 
three tighter than before. 

In this letter, we take advantage of these new measurements 
of the unresolved fraction of the CXRB to put firm upper limits 
on the global accretion history of massive black holes at z > 5. 
The aim of our approach is not to exclude a particular model 
but rather to highlight the existence of some tension between the 
need for efficient and rapid accretion required by the observation 
of supermassive black holes already in place at z — 7 and the 
strict upper limit on the accreted mass of the whole high-z MBH 
population imposed by the very tiny CXRB unresolved fraction. 

We adopt a (h, Q m , Q A ) = (0.7, 0.3, 0.7) cosmology. 



2. Soltan's argument and the CXRB 

The so-called "Soltan's argument" (Soltan 1982) translates the 
observed radiation emission of AGNs integrated over the cosmic 
history of the Universe into mass accreted onto the putative su- 
permassive black-hole population. While the argument is usually 
expressed in terms of AGN luminosity functions, it is straight- 
forward to apply it in the context of background radiation in a 
given observed band, which can easily be done as follows. 



1 



Ruben Salvaterra et al.: Limits on the high redshift growth of massive black holes 



Let us assume that the average AGN X-ray spectrum has a 
power-law form L v oc E~ a , so that the comoving specific emis- 
sivity vs. redshift can be factorized as 



(1) 



j(E,z) = U\ — | f(z). 



Here /(z) is a function describing the redshift evolution of the 
emissivity. Neglecting absorption in the IGM, the contribution 
to the background at observed energy Eo due to sources located 
at redshifts z > z is 
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where E - £\)(1 + z). If z » ZmA - 0.33 (the matter-A equiva- 
lence redshift), we can neglect the cosmological constant energy 
density in the line element, and by assuming f(z) = (1 + zY y 
derive the normalization emissivity 



J* 



* = 47r7 £o ^Q,f (Eo/E+na + y + 3/2)(l + i)^+ 3 / 2 . 



(3) 



Now let us consider the standard Soltan's argument. The co- 
moving mass density accreted onto MBHs within a given I is 



Paccfe) = 5— 
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dEj(E,z), 



(4) 



where e is the mass-radiation conversion efficiency. It is worth 
noting that at high redshifts p acc could well be significantly lower 
than the total mass density locked in MBHs. We now relate the 
bolometric emissivity to the emissivity in a given energy band, 
[E m - Em], by introducing a bolometric correction fx 



dEj(E,z)j J o dEj(E,z). 



(5) 



By substituting eq. 3 into eq. 4 and integrating, we finally obtain 
(1-e) E J Eo (a + r + 3/2) 



Paccfe) = 47T- 



ec 3 fx(l-a) (y + 3/2) 
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which is valid for a + 1 . The above formula allows us to esti- 
mate the maximum mass accreted onto MBHs within any given 
redshift interval, that contributes at any specified level to the ob- 
served background. 

3. Results 

We are interested in placing a firm upper limit on the mass 
accreted onto MBHs prior to z > 5 by considering the unre- 
solved fraction of the CXRB. The unresolved CXRB is well- 
described by a single power-law with a very hard photon-index 
(=* 0.1 + 0.7) and a flux of 2.5+} * x 1(T 12 erg s _1 cnr 2 deg" 2 
in the 1.5-7 keV energy band. Comparing the measured un- 
resolved CXRB to the AGN population model by Gilli et al. 
(2007), Moretti et al. (2012) find that most of the flux at =* 1.5 
keV can be accounted for by faint, z < 5 sources, but that 
their model falls short for E > 3 keV, suggesting that there is 
a larger population of Compton thick sources at moderate red- 
shifts (z =2 2, see fig. 10 in Moretti et al. 2012). The maximum 
allowed (1-er error) intensity of the unresolved CXRB at 1 .5 ke V 
is £1.5/1 5 0.47 X 10 -12 erg s _1 cirT 2 deg~ 2 . This value reduces 
to 0.21 x 10~ 12 erg s _1 cirT 2 deg~ 2 when the contribution of the 
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z < 5 faint sources modeled by Gilli et al. (2007) is taken into ac- 
count. Eq. 6 derived in the previous section can be used to readily 
translate these limits into constraints on the total accreted mass 
density. To model the average Type I-Type II composite AGN 
spectrum, we adopt the spectral energy distribution (SED) pro- 
posed by Sazonov, Ostriker & Sunyaev (2004). In the 2-10 keV 
range, the spectrum is well-approximated by a power-law with 
a =* 0.25, with a bolometric correction fx ~ 0.04. Assuming 
e — 0.1, for£o = 1.5 keV eq. 6 give£] 

feW <3.4xI0«(^)(f)(i^) 02I M B Mpc- 

In general, our mass density limit is valid for any X-ray emitting 
population that is not resolved above z in the deepest X-ray sur- 
veys. For our purposes, the above limit is valid in particular for 
z > 5 as no AGNs above such a redshift limit has been identified 
in the 4 Msec CDF-S (Xue et al. 201 1). 

In fig. 1, we plot the limits obtained by adopting the Moretti 
et al. (2012) results described above. The dark shaded region is 
excluded since it would imply an unresolved CXRB at 1.5 keV 
that is higher than observed. The more stringent light-shaded 
region is obtained by subtracting the contribution of faint or 
absorbed sources or both at z < 5. At z = 5, the maximal 
possible accreted mass along the cosmic evolution of MBHs is 
Pace < 1.4 x 10 4 M Mpc- 3 (p acc < 0.66 x 10 4 M Q Mpc~ 3 sub- 
tracting the faint sources at lower redshifts). 

The limits obtained above provide strong constraints on the 
models of the formation and evolution of MBHs in the early 
Universe. The accreted mass density is one of the most direct 
predictions of semi-analytical models (see Volonteri 2010 for 
a review) and can be directly compared with our constraints. 
We compare our limit on p acc to simple models (Volonteri & 
Begelman 2010), which assume that all MBHs shine at a fixed 
fraction, 30%, of the Eddington luminosity and that accretion ac- 
tivity is major-merger driven and self-regulated by the host, as- 
suming an unvarying relation to the velocity dispersion at all red- 
shifts. Gravitational wave recoil is neglected. In spite of the sim- 
ple description of the physical processes leading to the growth 
of the seeds (e.g., assumptions of a constant Eddington ratio 
for all objects at all redshifts and of the local scaling relations 
being established in the high redshift Universe), these models 
are rather successful in reproducing the observed AGN popula- 
tion. As an example, they reproduce the local MBH mass den- 
sity, and the observed AGN bolometric luminosity function at 
low-to-intermediate redshifts. We stress that the accreted mass 
density should be considered a qualitative result, as Volonteri 
& Begelman (2010) did not attempt to model in detail the evolu- 
tion of the MBH population (e.g., they assumed a fixed accretion 
rate for all MBHs at all times, which is clearly an oversimpli- 
fication), but instead whished to estimate whether a particular 
class of seeds produced a physically reasonable MBH popula- 
tion. Despite the qualitative nature of the theoretical constraints, 
these models are very useful in putting the strict limit on the ac- 
creted mass density implied by the CXRB in the context of MBH 
evolution studies. 

Models involving stellar-size Population III seeds may have 
difficulties in accounting for the existence of a MBH as mas- 
sive as M > 10 9 M at z = 7.1, as observed by Mortlock et 
al. (2011) (see discussion in Petri, Ferrara, Salvaterra 2012). A 
model involving stellar-sized seeds is below our limits (lower 
long-dashed curve in Fig. 1), but fixing the accretion rate at 30% 



1 We note that (a + y + 3/2)/(y + 3/2) ^ 1 for reasonable values of 
the cosmic evolution y. 



2 



Ruben Salvaterra et al.: Limits on the high redshift growth of massive black holes 



10 s F 



I 

o 
a 

2 

o 
2 



7) 

c 

CD 

■a 

w 
w 
to 



10" r 



1000 - 



XI 
<u 
-<-> 
<u 
u 
o 
a 

CO 



100 




7 8 
redshift 



Fig. 1. Limits on the density of accreted mass onto MBHs at 
z > 5 derived from the unresolved fraction of the CXRB ob- 
served at 1 .5 keV. Dark shaded area refers to the maximum al- 
lowed CXRB intensity and light shaded area is the limit once 
the contribution of lower-z AGNs is taken into account. For il- 
lustration, four curves show the accreted mass density of models 
of formation and evolution of MBHs presented in Volonteri & 
Begelman (2010). These mass densities should be considered 
qualitative, rather than quantitative, estimates, but they provide 
the typical range found when assuming a fixed accretion rate for 
all MBHs and self-regulated growth (see also Treister et al. 201 1 
for additional examples). The dotted curve refers to a model that 
induces an early reionization, and the short-dashed curve to a 
model that barely reproduces the mass function of z ~ 6 quasars. 
These two curves assume massive MBH seeds and /ndd = 0.3. 
The dot-dashed curve is analogous to the short-dashed model, 
but assumes the distribution of /em given by Merloni & Heinz 
(2008). The long-dashed curve is based on Population III rem- 
nants, and, while being consistent with the CXRB constraint, is 
unsuccessful in assembling 10 9 M MBHs by z — 6. 



of the Eddington luminosity, it is unable to explain the presence 
of the population of z — 6 quasars (Willott et al. 2010). If we 
increased the accretion rate to 100% Eddington, to account for 
z = 6 quasars, this model would overproduce the total unre- 
solved CXRB (see Treister et al. 2011 for an example of this 
case). 

Massive seeds, e.g., the so-called "quasistars" (see Volonteri 
2010 and references therein), seem to be a more viable option to 
explain the observed population of high-z quasars. This class of 
models accounts for a population of MBHs with ~ 10 9 M at 
z - 6-7. Volonteri & Begelman (2010) provide two observa- 
tionally limited cases for the efficiency of the formation of mas- 
sive seeds. The results in terms of p acc are also shown in fig. 1. 
The high efficiency model (upper dotted curve) results in a very 
early reionization of the IGM, while the low efficiency one (mid- 
dle short-dashed curve) barely succeeds in assembling enough 
MBHs as massive as 10 9 M by z — 6. Direct collapse models 
with fixed accretion rates all exceed our limits. The model with 



the higher efficiency would overproduce the total unresolved 
CXRB at z - 7.2, while the low efficiency model is above the 
more stringent limit based on the subtraction of lower-resdhift 
sources, for z > 6.3. Similar results are found by considering the 
MBH growth models of Volonteri, Lodato & Natarajan (2008), 
and Agarwal et al. (2012), among others. 

We stress that our result does not directly favour one seed 
model over another, but highlights the strong constraints on the 
average accretion rate of the high-z MBH population, as a whole. 
Assuming an initial mass density (po), and a mean Eddington 
raticQ(/EddX the accreted mass density is 



PaccM = P0 



exp / E dd 

t e 



1 
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where t - cr T c/(4nG m p ) = 0.44 Gyrs (c is the speed of light, 
(jj is the Thomson cross-section, and m p is the proton rest 
mass). This approximation requires that all seeds have similar 
masses and form roughly at the same time. Additionally, such a 
generic argument does not take into account any self-regulation 
or feedback effect that limits the MBH growth. Keeping these 
caveats in mind, for plausible values of po — 10 - 1000 M Q 
Mpc~ 3 , the average /ndd must be less than 0.1-0.3 at z > 5. 
This is in line with lower-redshift results that the distribution of 
Eddington rates of z — 2 - 4 luminous quasars is dominated by 
sub-Eddington sources (Kelly et al. 2010). With such an average 
/Edd, the M = 2 x 10 9 M black-hole observed by Mortlock et 
al. (201 1) at z = 7.1 would require a seed black-hole of mass ex- 
ceeding 2x 10 7 M . We conclude that the most massive MBHs at 
very high redshift cannot accrete at the average Eddington ratio 
(see also Trakhtenbrot et al. 201 1). Models require rates close to 
Eddington to explain the high-mass end of the mass function of 
quasar-powering MBHs (Natarajan & Volonteri 201 1). 



4. Possible solutions 

Massive black-hole growth models seem to predict p acc above 
our observational limit. One may argue that the bulk of the 
accreted mass is missing in the small volume sampled by the 
CDFs. A simple estimate implies that objects rarer than ~ 10~ 6 
Mpc~ 3 are not present in the field. However, even assuming that 
all of these BHs were accreted up to ~ 10 8 M Q , the resulting 
Pace would be much lower than our limit, further increasing the 
disagreement with evolutionary models. 

On the theoretical side, models tested here are based on sim- 
ple assumptions. As an example, p acc is expected to be reduced 
by gravitational wave recoil, which we however neglect. We did 
check that in the considered models the effect on p acc is at most 
20%, hence that the inclusion of the recoil does not repre- 
sent a viable solution. A stronger impact, but still insufficient, 
is obtained by considering a non-fixed accretion rate. A model 
starting from massive seeds coupled to the empirical distribution 
of Eddington ratios derived by Merloni & Heinz (2008), shown 
in fig. 1 as a dot-dashed line, exceeds our limit at z < 6. 

The most promising solution relies on the possibility that the 
most massive black holes are able to maintain a high /gdd dur- 
ing their cosmic history, while lighter ones accrete at a much 
lower rate. Therefore, accretion must strongly depend on either 
the MBH mass (or most likely the host mass as in 'selective ac- 
cretion', Volonteri & Stark 2011) or environment (di Matteo et 
al. 2012), or must have a low-mass cut-off ("global warming", 



2 Which can be considered a combination of accretion rate and duty 
cycle, see, e.g., Tanaka, Perna, & Haiman 2012. 



3 



Ruben Salvaterra et al.: Limits on the high redshift growth of massive black holes 



Tanaka et al. 2012). This is in line with the observational evi- 
dence that high redshift quasars seem to be powered by MBHs 
that are 'over-massive' for a fixed galaxy property with respect 
to their counterparts at z — (e.g., Wang et al. 2010). Willott et 
al. (2010) instead find that either many massive galaxies at z — 6 
do not have MBHs, or that these MBHs are less massive than 
expected assuming that MBHs are roughly 1/1000 of the host 
stellar mass. This suggests overall that while some (most-likely 
the most massive) MBHs can grow above today's correlations, 
most of them should be less massive than expected from local 
relations (cf. Volonteri & Stark 201 1). 

5. Conclusions 

We have placed firm upper limits on the global accretion his- 
tory of MBH at z ^ 5 by taking advantage of the measurement 
of the unresolved fraction of the CXRB provided by Moretti et 
al. (2012). The maximum allowed unresolved CXRB intensity 
observed at 1.5 keV implies a maximum accreted-mass density 
onto MBH at z > 5 of p acc < 1.4 x 10 4 M Mpc~ 3 . Considering 
the contribution of lower-z AGN^ (Gilli et al. 2007), this limit 
reduces to p acc ~ 0.66 x 10 4 M Q Mpc~ 3 . This value translates 
into < 1 HI ionizing photon per baryon produced by accretion 
onto MBHs at z — 6, confirming the common wisdom that hy- 
drogen reionization is driven by stellar-like sources (cf. Haardt 
& Madau 2012). 

It is important to stress that this calculation is a strict up- 
per limit to the accreted mass onto MBHs at z > 5, as most of 
the unresolved CXRB could well be attributed to faint sources 
at lower redshift that are not accounted for by the Gilli et al. 
model. That the spectrum of the unresolved CXRB determined 
by Moretti et al. (2012) is very hard, alone seems to exclude 
there be a significant contribution from AGNs at z > 5. For such 
a population, the Compton reflection roll-over (rest-frame en- 
ergy ^ 30 keV) would fall within the observed energy band, 
resulting in an emission much softer than the observed unre- 
solved CXRB. Therefore, a significant contribution from high-z 
sources would result in an even harder spectrum of the still un- 
accounted fraction of the CXRB, possibly at odds with current 
population-synthesis models. 

A possibly stricter upper limit on p acc at z - 6.5 can be 
obtained by the stacking analysis of the X-ray emission of i- 
dropouts selected by Bouwens et al. (2006) in the CDF-S. In 
contrast to Treister et al. (2011), Willott (2011), Fiore et al. 
(2012), and Cowie et al. (2012) did not find any evidence of 
X-ray emission. The flux limit in the observed hard X-ray band 
derived from the stacking analysis corresponds to a MBH mass 
density of p acc < 0.4 x 10 4 M Mpc~ 3 . A much stronger, about 
ten times lower, upper limit is obtained from the flux limits ob- 
tained using the more sensitive soft X-ray band, though here 
absorption can play a decisive role. These limits are tighter than 
ours at the same redshift, although we believe that our results are 
less subject to biases and assumptions than those derived from 
the stacking analysis. In particular: i) the stacking analysis relies 
on the corrections for incompleteness, photometric redshift mea- 
surements, and dust absorption of the Bouwens et al. sample; ii) 
the flux limit of the /-dropouts implicitly introduces a lower limit 
on the MBH mass probed (that we estimate to be a few times 10 6 
M ); and iii) the stacking analysis results strictly refer only to the 
AGN activity in a narrow redshift range. The limits we place on 

3 Adopting the speculative model of Moretti et al. (2012) which is 
able to explain the very hard shape of the unresolved 1.5-7 keV CXRB, 
the limit on p acc further decreases by a factor of ^ 2. 



p acc are instead unaffected by any of these effects since the back- 
ground intensity directly measures the time-integrated accreted 
mass. 

We investigate how much our assumptions on the emission 
properties of high-z MBHs affect our results. The adoption of 
a much flatter SED, e.g. with a = 0.9, will result in limits that 
are 1.8 less stringent than those shown in fig. 1. On the other 
hand, as shown by Marconi et al. (2004), the fraction of light 
emitted by AGNs in the rest-frame 2-10 keV band increases 
with decreasing bolometric luminosities, so that fx = 0.04 can 
be considered a conservative choice. We note that at z > 5 the 
observed 1 .5 keV photons has been emitted at energy > 9 keV, 
where absorption can not be very strong even for high intrinsic 
column-densities. Thus, even assuming that high-z AGNs are 
heavily obscured, our limits cannot increase by much. Finally, 
our adopted value of the accretion efficiency, e = 0.1, is in- 
between the range of allowed values. For Schwartzschild BHs 
(e = 0.057), the limits would be 1.7 higher, while for maximal 
rotating BHs (e = 0.42) they are a factor of six lower. 

The bottom line of our analysis is that there is some tension 
between the need for efficient and rapid accretion required by the 
observation of SMBHs already in place at z = 7.1 (Mortlock et 
al. 201 1), and the strict upper limit on the accreted mass derived 
from the CXRB. Therefore, accretion must be very efficient for 
the most massive BHs during their lifetime, but sub-Eddington 
for most of the AGN population. 
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